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Using a sample of 31.3 x 10 6 BB pairs collected with the Belle detector at the T(4S) resonance, we 
make the first observation of the charged B meson decay to XcO and a charged kaon. The measured 
branching fraction is B(B + — > \coK + ) = (6.0^ 2 'g ± 1.1) x 1CP 4 , where the first error is statistical, 
and the second is systematic. 
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Two-body decays of B mesons with a charmonium par- 
ticle in the final state have recently received substantial 
attention due to their sensitivity to CP violation in B 
system. The production rate of charmonium states in 
b — ► ccs transitions also provides valuable insight into 
the dynamics of strong interactions in heavy meson sys- 
tems. For instance, although the production of the XcO 
P-wave ++ state in B decays vanishes in the factoriza- 
tion approximation as a consequence of spin-parity and 
vector current conservation, it is possible if there is an 
exchange of an additional soft gluon pU|. At present, 
only upper limits on the B — > XcoK branching fractions 
exist m. In this Letter, we report the first observation of 
the B + — > XcoK + decay mode. The analysis is performed 
using data collected with the Belle detector at the KEKB 
asymmetric energy e + e~ collider jjj. The data sample 
consists of 29.1 fb _1 taken at the T(45) resonance con- 
taining 31.3 x 10 6 BB pairs, and 2.3 ftr 1 taken 60 MeV 
below the T(45) resonance to perform systematic studies 
of the e + e" — * qq background. 

The Belle detector || is a large-solid-angle magnetic 
spectrometer that consists of a three-layer silicon ver- 
tex detector (SVD), a 50-layer central drift chamber 
(CDC) for charged particle tracking and specific ioniza- 
tion measurement (dE/dx), an array of aerogel thresh- 
old Cerenkov counters (ACC), time-of- flight scintillation 
counters (TOF), and an array of 8736 CsI(Tl) crystals 
for electromagnetic calorimetry (ECL) located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux return located outside the coil is 
instrumented to detect Kl mesons and to identify muons 
(KLM). Electron identification is based on a combina- 
tion of CDC dE/dx measurements, the response of the 
ACC, and the position, shape and energy deposition of 
the associated ECL shower. We use a Monte Carlo (MC) 
simulation to model the response of the detector and de- 
termine acceptance ||. 

The analysis of the B + — + \ c qK + decay is performed in 
the framework of a general study of the B + — > K + h + h~ 
decay, where h stands for either pion or kaon H (charge 
conjugation is implied throughout this Letter). We re- 
construct the XcO meson in the decay modes XcO — ► 7r + 7r~ 
and XcO - * K + K~ . Charged tracks are required to sat- 
isfy requirements based on the average hit residual and on 
their impact parameters relative to the interaction point. 
We require that the transverse momentum of the track 
be greater than 0.1 GeV/c to reduce low momentum 
combinatorial background. For charged hadron identi- 
fication, we use a combination of CDC dE/dx measure- 
ments, flight time measured in TOF, and the response 
of the ACC. We select kaon candidate tracks with a set 
of criteria that has about 90% efficiency, a charged pion 
misidentification probability of about 8%, and a negligi- 
ble contamination from protons. We reject tracks that 
are positively identified as electrons. 

We reconstruct B mesons by combining a XcO with 



a charged kaon. The candidate events are identi- 
fied by their center-of-mass (cm.) energy difference, 
AE = (J^i Ei) — Et, and the beam constrained mass, 
Mb c = \jEl~ {J2iPi) 2 i where E^ = y/s/2 is the beam 
energy in the cm. frame, and pi and Ei are the 
cm. three-momenta and energies of the candidate B 
meson decay products. We select events with M& c > 
5.20 GeV/c 2 and \AE\ < 0.2 GeV, and define a signal 
region of \M bc - M B \ < 9 MeV/c 2 and \ AE\ < 0.04 GeV 
and two AE sideband regions —0.08 GeV < AE < 
-0.05 GeV and 0.05 GeV < AE < 0.15 GeV g. 

To suppress the large combinatorial background, which 
is dominated by the two-jet-like e + e~ — > qq contin- 
uum process, we use variables that characterize the event 
topology. We require |cos6*thr| < 0.80, where # t hr is the 
angle between the thrust axis of the B candidate and that 
of the rest of the event. This eliminates 83% of the con- 
tinuum background and retains 79% of the signal events. 
We also define a Fisher discriminant, T , that includes the 
production angle of the B candidate, the angle of the B 
candidate thrust axis with respect to the beam axis, and 
nine parameters that characterize the momentum flow in 
the event relative to the B candidate thrust axis in the 
cm. frame ||. We impose a requirement on T that re- 
jects 79% of the remaining continuum background and 
retains 74% of the signal. In the case where the XcO is 
reconstructed in the K + K~ mode, the continuum back- 
ground is much smaller and a looser requirement that re- 
jects 53% of the continuum background with about 89% 
efficiency for the signal is used. 

As shown in Ref. JjJ, the background to the K + ir + -K~ 
final state from BB decays is dominated by decays of 
the type B + — > \K + ir~]n + , where [-ftT + 7r~] denotes an 
intermediate state that can decay into the K + ir~ such 
as K*°(892) or D°. To suppress this type of background, 
we require that the invariant mass of the K + tt~ system 
be greater than 2.0 GeV/c 2 . For the K+K+K' final 
state, we require that the invariant mass for both K + K~ 
combinations be greater than 2.0 GeV/c 2 to suppress the 
background from charmless B decays. 

We select all K + h + -k~ (K + K + K~) combinations 
from the B signal region that satisfy the selection crite- 
ria described above and have a 7r + 7r~ (K + K~) invariant 
mass in the range 3.2 GeV/c 2 < M(h+h~) < 3.8GeV/c 2 . 
The resulting tt + tt~ and K + K~ invariant mass spectra 
are shown in Fig. [l| Since in the case of the K + K + K~ 
final state there are two kaons with the same charge, 
we distinguish the K + K~ combinations with the smaller 
and larger invariant masses. Here, we only plot the 
larger of the two possible combinations and require the 
invariant mass of the smaller combination to be less than 
3.35 GeV/c 2 . The peaks near 3.4 GeV/c 2 in Figs. |(a) 
and |l|(b) are identified as the XcO meson. The peak at 
3.69 GeV/c 2 in Fig. |l|(a) corresponds to the '0(25) meson 
from the decay mode B+ -> tJj(2S)K + , 0(25) -> (J, + fi~ 
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FIG. 1. The (a) tt + -k~ and (b) K + K~ invariant mass spec- 
tra. Open histograms correspond to events from the B signal 
region, and hatched histograms correspond to events from the 
AE sidebands. The curves are fits to the signal data. 



with muons misidentified as pions. The hatched his- 
tograms in Fig. [l] correspond to the events from the AE 
sidebands plotted with a weight of 0.62. 

These 7r + 7r~ and K + K~ spectra are fitted in the range 
3.20 GeV/c 2 < M{h+h-) < 3.62 GeV/c 2 to the sum of a 
constant for the background and a Breit-Wigner function 
convolved with a Gaussian resolution function for the 
signal. The width of the resolution function is fixed at 
11.0 MeV/c 2 as determined from a fit to the J/ip peak 
in the invariant mass spectrum. The full width of 

the Breit-Wigner function is fixed at the world average 
Xco width of 14.9 MeV @. The results of the fit are 
given in Table |. The statistical significance, S, of the 
signal, in terms of the number of standard deviations, 
is calculated as \f— 2 ln(£o / £ max ) > wnere £max and Cq 
denote the maximum likelihood with the nominal signal 
yield and with the signal yield fixed at zero, respectively. 

The peak position in the K + K~ spectrum is found to 
be shifted below the PDG value f| by 25 ± 10 MeV/c 2 : 
see Fig. [j](b) and Table |. Although this shift is consistent 
with a statistical fluctuation, we note that, as shown be- 
low, we observe evidence for a non-resonant-like signal in 
the B+ ->■ K+K+K- final state. As a result, the K+R- 
invariant mass distribution in the Xco region could be 
distorted by the effects of interference with an amplitude 
not related to the B + — > \cqK + . Because of this uncer- 
tainty, we base our branching fraction measurement on 
the \c0 — > 7r + 7r~ decay mode only. 

For the branching fraction calculation, we normalize 
our results to the observed B + — ► J/-0A+, J/ip — > 
signal. This removes systematic effects in the particle 
identification efficiency, charged track reconstruction ef- 
ficiency and the systematic uncertainty due to the cuts 
on event shape variables. To avoid additional systematic 
uncertainty in the muon identification efficiency, we do 



TABLE I. Results of the fit to the tt+tT 
H + invariant mass spectra. 
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not use muon identification information for J/tp recon- 
struction. Instead, we apply the same pion-kaon sepa- 
ration requirement for muons from the J /ip as for pions 
from the Xco- The feed-across from the J/ip — > e + e~ 
submode is found to be negligible (less than 0.5%) after 
the application of the electron veto requirement. 

To determine the B + — > XcoK branching fraction, 
we use the signal yield obtained from the fit to the tt + 7t~ 
invariant mass spectrum. The number of J /%pK + signal 
events is determined from the fit to the invariant 
mass spectrum (see Table |) . Combining all the relevant 
numbers from Table § and using the intermediate branch- 
ing fractions of B(xcO — > 7r + 7r~) = (5.0 ±0.7) x 10~ 3 and 
B(J/4) -► fM+fi-) = (5.88 ±0.10) x 10- 2 0, we find the 
ratio of branching fractions: 



B(B+ -> XcoK + ) 
B(B+ -» J/i>K+) 



0.60 J 



0.18 



±0.05 ±0.08, 



where the first error is statistical, the second is sys- 
tematic, and the third is due to the uncertainty in the 
XcO — > 7r + 7r~ branching fraction. Here the statistical er- 
ror includes the errors on the number of signal XcoK + 
and J/ipK + events. The systematic error consists of 
the uncertainty in the J/ip — > (J- /jT branching frac- 
tion (1.7%) and the uncertainty in the background pa- 
rameterization in the fit to the it + tt~ spectra (7.8%). 
Using the world average value of B(B + — > J/ipK + ) = 
(10.0 ± 1.0) x 10~ 4 jl(|, we translate our measurement 
into the branching fraction: 



B(B+ -» XcoK + ) 



(6.0±f ± 1.1) x 10~ 4 , 



where the first error is statistical, and the second is 
the total systematic error including uncertainties in the 
XcO —> 7r + 7r~ and B + — » J /%jjK + branching fractions. 

Figure || shows the projections of the AE signal 
bands for the selected B + — > XcqK + candidates with 
\M(h+h-) - M Xc0 \ < 0.05 GeV/c 2 and 5.270 GeV/c 2 < 
M bc < 5.288 GeV/c 2 . The hatched histograms in 
Fig. U correspond to the events from the Xco mass side- 
bands defined as 0.07 GeV/c 2 < \M(h + h~) - M Xa0 \ < 
0.17 GeV/c 2 and plotted with a weight of 0.50. For the 
A' + 7r + 7r~ final state, the distribution for the sideband 
events is consistent with background. In contrast, in 
the three charged kaon final state we observe a substan- 
tial signal for events in the Xco mass sidebands. This 
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FIG. 2. The AE distributions for B + -> XcoK + candi- 
dates; (a)-XcO - * 7T + 7r _ and (b)-XcO — ► K + K~ . The hatched 
histograms correspond to the XcO mass sidebands. The solid 
lines display the signal plus background combined shape. The 
dashed lines correspond to the background shape only. 

is the evidence for non-resonant-like B + — » K + K + K~ 
decays that may be responsible for the shift in the 
Xco K + K~ mass peak discussed above. A study 
of the full K + K + K~ Dalitz plot also supports this 
conclusion. We represent the background shape in AE 
with a linear function and restrict the fit to the range 
-0.1 GeV < AE < 0.2 GeV {§. The signal shape 
is parameterized by a sum of two Gaussians with the 
same mean. The AE shape for the signal is determined 
from the B + — ► D°ir + events. The B + — > XcoK + where 
Xco -> 7T+7T- (if+if-) signal yield of 16.4±p (9.7^'j) 
events obtained from the fit to the AE distributions 
agrees with that obtained from the fit to ir + ir~ (K + K~) 
invariant mass spectrum. 

To cross-check the result, we also reconstruct the Xco 
meson in the K + K~ir + n~ final state. We reduce the 
large combinatorial background by using events where 
at least one Kn pair has an invariant mass within 
75 MeV/c 2 of Mk*- We also apply the tighter require- 
ment |cos#thr| < 0.6 to suppress the continuum back- 
ground. Figure || presents the K* (892)Kir invariant 
mass spectrum for the selected events from the B sig- 
nal region shown by open histogram and for events from 
the AE sidebands shown by hatched histogram plot- 
ted with a weight of 0.62. For the branching fraction 
calculation we normalize our result to the signal ob- 
served in the B+ -► J/iPK+, J/xjj -> K*{892)K-ir+ 
decay mode. Figure [| shows the AE distributions for 
B + — > XcqK + and B + — ► J/ipK + candidates. From 
the fit to the AE distributions we find 21.9 ± 6.3 and 
18.0 ± 5.2 signal events, respectively. The statistical sig- 
nificance of the B + — ► XcoK + signal in this mode is 4.3<r. 
Using the intermediate branching fractions of B(y c o — ► 
K*°(892)K-ir+ + c.c.) = (1.2 ± 0.4) x 10" 2 |Mj and 
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FIG. 3. The K* (892)Kn invariant mass spectrum for 
events from the B signal region (open histogram) and for 
events from the A_B sidebands (hatched histogram). 
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FIG. 4. The AE distributions for (a) B + -» XcoK + , 
Xco -> K*° (892) K~tt+ candidates and (b) B+ -*• J/ipK+, 
J/ip — » K*° (892)K ~7r + candidates. The hatched histogram 
in (a) corresponds to events from the XcO mass sidebands. 
The solid lines display the signal plus background combined 
shape. The dashed lines correspond to the background shape 
only. 



B(J/i/j -> K*° (892) R-TT+ + c.c.) = (6.3±2.1) x 10~ 3 @, 
we determine the ratio of branching fractions: B(B + — > 
XcoK+)/B(B+ -> J/iPK+) = 0.64 ±0.26 ±0.30, where 
the first error is statistical, and the second consists of the 
uncertainty in the secondary branching fractions. The 
obtained number is in agreement with that determined 
from the Xco 7r + 7r - decay mode. We do not include 
this result in the final value for the B + — > XcoK + branch- 
ing fraction because of the large systematic uncertainty. 

In summary, we report the first observation of the 
B + — * XcqK + . The statistical significance of the sig- 
nal is 6cr when the Xco ~~ * 7r + 7r _ and Xco ~~ * K*Kit 
modes are combined. The measured branching fraction is 
B{B+ -> XcqK+) = (6.0t 2 ;g±l.l)xl0~ 4 which is compa- 
rable to those for the B + — » J/tpK + and B + — » XciK + 
decays. This provides evidence for a significant nonfac- 
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torizable contribution in B to charmonium decay pro- 
cesses. The result reported here supersedes the previous 
value based on 21.3 fb _1 as reported in Ref. flT3| . 
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this reference. From Fig. 9 therein, we estimate 120 ± 13 
signal events in the 7^*°(892)_K' _ 7r + final state. Assuming 
the same reconstruction efficiency as for the K + K~tt + tv~ 
final state, we obtain B(J/ip -+ K*° (892)K~ ty + + c.c.) = 
(6.3 ±2.1) x 10 -3 . 
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